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day-night (D-N) effect for the Super-Kamiokande detector in the case of the MSW 

v e — ► VuM transition solution of the solar neutrino problem. The one-year averaged 
probability of survival of the solar v e crossing the Earth mantle, the core, the inner 2/3 
i-£h \ of the core, and the (core + mantle) is calculated with high precision (better than 1%) 

using the elliptical orbit approximation (E0A) to describe the Earth motion around 
the Sun. Results for the survival probability in the indicated cases are obtained for 
a large set of values of the MSW transition parameters Am 2 and sin 2 29y from the 
"conservative" regions of the MSW solution, derived by taking into account possible 
relatively large uncertainties in the values of the 8 B and 7 Be neutrino fluxes. Our 
results show that the one- year averaged D-N asymmetry in the v e survival probability 
for neutrinos crossing the Earth core can be, in the case of sin 2 29y < 0.13, larger than 
the asymmetry in the probability for (only mantle crossing + core crossing) neutrinos 
by a factor of up to six. The enhancement is larger in the case of neutrinos crossing 
the inner 2/3 of the core. This indicates that the Super-Kamiokande experiment 
might be able to test the sin 2 29y < 0.01 region of the MSW solution of the solar 
neutrino problem by performing selective D-N asymmetry measurements. 



1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria. 



1 Introduction 

The Earth effect in neutrino propagation is a direct consequence of neutrino oscillations in 
matter |l|. If the neutrino mass spectrum is not degenerate and if neutrino mixing takes 
place in vacuum, the neutrino propagation is sensitive to the matter distribution along the 
propagation path and the probability to detect at Earth a v e produced in the Sun is a 
function of the detection time because at night the Sun is below the horizon and the solar 
neutrinos cross the Earth reaching the detector. 

The observation of the Earth, or day- night (D-N), effect would be a proof of the validity 
of the MSW solution of the solar neutrino problem. In its simplest version the MSW 
mechanism involves matter - enhanced two-neutrino transitions of the solar v e into an 
active neutrino v^ or u T , u e —>■ i/j t ), while the v e propagates from the central part to the 
surface of the Sun. The v e —>■ v^tr) transition probability depends in this case on two 
parameters: Am 2 and sin 2 26V, where Am 2 > is neutrino mass squared difference and 9 V 
is the angle characterizing the neutrino mixing in vacuum. When the v^m produced in the 
Sun due to the MSW transitions and the solar u e not converted into u^ T ) in the Sun cross 
the Earth, further v e — > v^m transitions and/or v e regeneration due to the inverse MSW 
process ^^m - > v e can take place. This can lead to a difference between the signals caused 
by the solar neutrinos in a solar neutrino detector during the day and during the night, i.e., 
to a D-N asymmetry in the signal. No other mechanism of depletion of the solar u e flux 
proposed so far can produce such an effect. 

The MSW solution of the solar neutrino problem and the D-N effect related to it have 
been extensively studied. The latest solar neutrino data can be described for values of the 
parameters Am 2 and sin 2 20y belonging to the intervals @ 



3.6 x 10~ 6 eV 2 < Am 2 <9.8 x KT 6 eV 2 , 

(1) 
4.5 x 10~ 3 < sin 2 28 v < 1.3 x 10" 2 , 

or 



5.7 x 1(T 6 eV 2 < Am 2 < 9.5 x 10" 5 eV 2 , 

(2) 
0.51 < sin 2 26 v < 0.92, 

assuming the solar v e MSW transitions are into active neutrinos. The intervals ([]]) and (^) 
correspond to the nonadiabatic (small mixing angle) and to the adiabatic (large mixing an- 
gle) solutions of the solar neutrino problem. They have been obtained using the predictions 
of the solar model of Bahcall and Pinsonneault from 1995 [[| for the fluxes of the different 
solar neutrino flux components: pp, pep, 7 Be, 8 B and CNO. This model takes into account 
heavy elements diffusion and was recently shown to be in very good agreement with the 
latest and more precise helio-seismological observations Q. 

The possible solution values of Am 2 and sin 2 26y are important in the calculations of 
the D-N effect because they determine the magnitude of the effect. A more conservative 
approach to the determination of the MSW solution regions of values of Am 2 and sin 2 29y 



allows for, e.g., the spread in the predictions for the 8 B, 7 Be, etc. fluxes in the contemporary 
solar models, the possible changes in the predictions associated with the uncertainties in 
the relevant nuclear reaction cross-sections which are used as input in the flux calculations, 
etc.Q. Such an approach || yields larger regions for both the nonadiabatic and the adiabatic 
MSW solutions 01: 



3.0 x 1(T 6 eV 2 < Am 2 < 1.2 x 10~ 5 eV 2 , 

(3) 
6.6 x 10~ 4 < sin 2 26 v < 1.5 x 1(T 2 , 

or 



7.0 x 10" 6 eV 2 < Am 2 < 1.6 x 10~ 4 eV 2 , 

(4) 

vr' n " 



0.30 < sin 2 26 v < 0.94. 



It is hardly conceivable at present that the values of Am 2 and sin 2 29y for which the MSW 
v e -^ v i>.(r) transition mechanism provides a solution of the solar neutrino problem can be 
very different from those shown in Fig. 1 and given in equations (|3|) and (f|). 

The D-N effect was widely studied g 0, g |, [TO], [TTj, |T2], [TJ, |TJ, |T§ . It was pointed out, 



in particular, that the effect can be significant for sin 2 29y > 0.01 for a large range of values 
Am 2 and sin 2 29y from the large mixing angle solution region, for which the Earth effect 
regenerates electron neutrinos increasing the event rate at night. At the same time results 



from Kamiokande II and III experiments [|Ty] did not reveal any significant Earth effect 



exceeding approximately 30%, thus excluding a rather large area of the region of values of 
Am 2 and sin 2 20y of the large mixing angle solution [13|. This area is already excluded by 



the existing mean event rate data from the Cl-Ar, Ga-Ge and Kamiokande experiments. 

The motivation for new quantitative theoretical studies of the D-N asymmetry comes 
from a number of considerations regarding the earlier studies and the experimental situation 
which is reasonable to expect in the near future. First of all the discovery of a D-N 
effect, even of a small one, should be an important confirmation of the neutrino oscillation 
hypothesis and will open a new chapter of physics, while the absence of any detectable 
effect should represent a strong constraint on the neutrino mass and mixing parameters 
relevant to the MSW solution. It is clear that in both cases quantitative predictions are 
required to assure a proper interpretation of experimental results. Second, most of the 
previous studies were primarily concerned with the nature of the effect || [| ||, [L0| . The 
relevant probability distributions and event rates were calculated for some representatives 
values of the parameters only, which is insufficient for a comprehensive understanding of the 
possible magnitude of the effect. Alternatively, in some of the articles on the subject the 
effect of the Earth on the oscillations of neutrinos was investigated for high energy neutrinos 
produced by cosmic rays or accelerators and the results obtained were extended to solar 



neutrinos ||, [17|, |TT], [L8| (see also, e.g., |T9|)- The possible application of the Earth effect to 



the problem of the reconstruction of the Earth internal structure (Earth Tomography) was 

2 Although the total fluxes of 8 B, 7 Be, etc. neutrinos depend on the physical conditions in the central 
part of the Sun, the spectra of the 8 B, CNO and the pp neutrinos are solar physics independent. 



also considered |[7], |2T|, ^] (see also || 0, [TTj] ) . More detailed results on the Earth effect 



have been obtained in (12j where the D-N asymmetry for the Super - Kamiokande detector, 
for instance, was calculated for 8 (4) sets of values of sin 2 29y and Am 2 from the region ([!]) 
((Q)) of the nonadiabatic (adiabatic) solution, as well as in [Q where the one year average 
iso - (D-N) asymmetry contours in the Am 2 - sin 2 29y plane for the Super - Kamiokande 
and SNO detectors have been derived. Third, the Kamiokande experiment, like the other 
solar neutrino experiments which have provided data so far, has a relatively low statistics. 
This problem will be overcome by new experiments such as Super - Kamiokande, SNO and 
ICARUS. They are expected to be much more sensitive to the D-N effect, allowing to exploit 
a wider neutrino parameter region where the Earth effect is not large, in particular, the 
region of the small mixing angle nonadiabatic solution. One of them, Super - Kamiokande, 
is already operating and is expected to produce soon new results. Finally, an improved 
statistics allows to develop more sophisticated tests enhancing the sensitivity to the effect. 

The present article represents the first of two articles aimed at providing comprehensive 
predictions for the D-N effect for the Super - Kamiokande detector in the case of the MSW 
v e -^ v h(t) transition solutions of the solar neutrino problem. These include: i) calculation 
of the Earth effect with a sufficiently high precision which can match the precision of the 
data on the effect to be provided by the Super - Kamiokande detector, ii) calculation of 
the magnitude of the effect for a sufficiently large and representative set of values of the 
parameters from the "conservative" regions (|3]) and (f|) of the MSW solution, iii) calculations 
of the one- year averaged D-N asymmetry for solar neutrinos crossing the Earth mantle, the 
core, the inner 2/3 of the core, and the mantle + core (full Earth) for the chosen large 
set of values of Am 2 and sin 2 29y with the purpose of illustrating the magnitude of the 
enhancement of the effect which can be achieved by an appropriate selection of the data 
sample, iv) calculation of the one year average deformations of the e~ spectrum by the 
Earth effect in the mantle, core, and in mantle + core for the selected representative set of 
values of Am 2 and sin 2 29y. 

The magnitude of the D-N effect is determined by the MSW probability of solar v e 
survival when the solar neutrinos cross the Earth to reach the detector, P®{v e — > u e ), as 
well as by the probability of solar u e survival in the Sun, -P©(f e — > z/ e ), describing the 
MSW effect on the solar neutrino flux detected during the day. Thus, the first step in 
the realization of our program of studies of the D-N effect consists of sufficiently precise 
calculation of the one year average probability P®(v e — > u e ) in the cases of interest, i.e., 
for neutrinos crossing the mantle, the core, 2/3 of the core and the mantle + core when 
traversing the Earth. In the present article this first step is accomplished. 

2 Calculating the Earth Effect 

Given that Super - Kamiokande detector should have an event rate of about 10 4 solar 
neutrino induced events per year (with a 5 MeV energy threshold), the expected statistical 
accuracy in the flux measurement is about 1%. This means that the numerical accuracy in 
computing the Earth effect has to be better than approximately 1% each time the effect is 
not smaller than 1%. 

Comprehensive and sufficiently accurate predictions require complex computations. In 
order to obtain detailed predictions for the magnitude of the D-N effect and how it changes 
with the change of the values of the two parameters Am 2 and sin 2 29y, we have chosen to 



calculate the effect for a large set of points from the MSW solution regions ([3]) and (f|). 
The values of Am 2 and sin 2 29y corresponding to these points are given in Table I. The 
points selected are distributed evenly in the "conservative" regions of the nonadiabatic and 
the adiabatic solutions (|3|) and (|J) 0. 

We have used the Stacey model |23| as a reference model for the Earth density distri- 
bution in our calculations. Obviously, this is not an unquestionable choice because new 
models were produced in the last years (see, e.g., p4|). However, the merit of the Stacey 



model when compared with alternative ones is its ability to describe the known details of 
the Earth radial density distribution, p(R), relevant for the Earth effect calculations. This 
is done through a set of polynomials which allows the use of very efficient numerical com- 
putation schemes. In the Stacey model the Earth is assumed to be spherical and isotropical 
with a radius of 6371 km. As in all Earth models, there are two main density structures 
in the Stacey model: mantle and core. The core has a radius of 3485.7 km and a mean 
density of approximately 11.5 gr/cm 3 . The mantle surrounding the core has a mean density 
of about 4 gr/cm 3 . Both of them have a number of density substructures. 

The Earth effect formalism was described by many authors and it was reviewed recently 
in ref. flI2|. The probability that an electron neutrino produced in the Sun will not be 



converted into v^t) when it propagates in the Sun and traverses the Earth on the way to 
the detector is given by [0] 



P e (^ e -> i/.) = P Q (u e -> i/.) + 1 2P °^p^ ( Pe2 _ S in 2 Ov), (5) 

COS Zuy 

where P<7){v e — > z/ e ) is the average probability of solar v e survival in the Sun and P e 2 is 
the probability of the z/ 2 — > u e transition after the u e have left the Sun. During the day, 
when the neutrinos do not cross the Earth, P e 2 = sin 2 6y and we have P®(v e — * v e) — 
P Q (z/ e — > v e ). For Earth crossing neutrinos at night P e 2 ^ sin 2 6y due to the MSW effect 
and -P®(z/ e — » Ve) 7^ P®{y e — > v e ). Throughout this study we use for pQ,(u e — > u e ) the 
expression derived in ref. |25j in the exponential approximation for the variation of the 
electron number density, N e , along the neutrino path in the Sun: 



P (z/ e -> v e ) = - + (- - P') cos 2# y cos 2d m , (6) 

where 8 m is the neutrino mixing angle in matter at the point of v e production in the Sun 
and 



-27rro^sin 2 20„ _ -2-ktq^- 

P' = ~ ~ (7) 

O^ Am 2 V ' 

is the analog of the Landau-Zener jump probability for exponentially varying N e , r being 
the scale height characterizing the change of N e along the neutrino trajectory. The proba- 
bility -P©(^ e — > v e ), eq. (0), was averaged over the region of 8 B neutrino production in the 
Sun. For values of sin 2 29y < 4 x 10~ 3 the probability Pq{v s — > u e ) was calculated following 



the prescriptions given in [2q . 



Let us discuss next the calculation of the probability P e2 . In contrast to P Q {v e — > v e ), the 
probability P e2 is not a constant in time because it is a function of the trajectory followed by 
the neutrinos crossing the Earth. The latter is determined by the instantaneous apparent 
position of the Sun in the sky. Given the fact that data are taken over a certain interval of 
time, the instantaneous P e2 in eq. (^) has to be replaced with its time averaged, < P e2 >, 
which requires the computation of the apparent solar trajectory in the detector sky. 

The exact solar trajectory is a complicated function of time which changes with the 
epoch of the year and the geographical location. Its computation is described in Appendix 
A, while for the present discussion it is enough to recall its main features. 

Looking at the Sun as at a star projected on a sphere centered on the detector location, 



its position is specified by a couple of spherical co-ordinates |27], |36|. Assuming the Earth 
is a spherical and isotropical body, the Sun's Nadir angle is enough to describe the Sun 
position in order to compute the probability P e2 . The Nadir angle h is the angle subtended 
by Sun's and Earth center directions as seen from the detector |, as depicted in Fig. 1. 
During the day h is greater than 90° (and P e2 = sin 2 #y), while at night h < 90°, the Sun 
rises or sets when h = 90° and a solar neutrino will cross the Earth center when h = 0°. 

The Sun's apparent trajectory in the sky with respect to the detector is a function 
of the day of the year d. Each night the Sun reaches a minimum Nadir angle, h m ^. In 
the course of the year /i m d has a maximum value, h m>max , and a minimum value, h m> min- 
Thus, given a point in the interior of the Earth seen (with respect to the detector) along 
a line of view of Nadir angle h, it will never be on the solar neutrino trajectory during 
the night if h < h mtmin ; it will be on a trajectory for some of the nights during the year if 
hm,min < h < h m ,max, while if h > h m ^ max it will be crossed by neutrinos every night. 

These limits and the apparent trajectory of the Sun, are functions of the detector loca- 
tion specified by its latitude, A^, and longitude. For a detector outside the tropical band, 
| Ad | > 23°27', h m ,max is reached at summer solstice, h mtmin is reached at winter solstice, 
and hm^min > 0° so that neutrinos received by such detectors never cross the Earth center. 
Instead, for detectors located inside the tropical band, h m ^ min = 0° and neutrinos cross- 
ing the Earth center are detected, while the minimum and maximum night Nadir angle 
are reached at epochs between equinoxes and solstices. At last for an equatorial detector 
h = 0° is reached just at equinoxes. 

It is well know that for a real time detector as Super - Kamiokande it is possible to 
compute h for the Sun at the time t of the detection of a given event, allowing the selection 
of events produced when the Sun is seen within a given Nadir angle interval [hi, /i 2 ]. The 
boundaries of this interval are crossed at times ti and t 2 . As will be shown further, this 
selection may be a feasible strategy to increase the sensitivity to the D-N effect. 

It follows from the above discussion that in order to take into account the Earth effect in 
the signal collected in the time interval [Ti, T 2 ] of the year by a given detector it is enough 
to replace P e2 in eq. fl5|) with its time averaged < P e2 > defined as: 



< Pei > = ^[ T2 dt5(h 1 < h(t) < h 2 )p e2 (h(t)), (8) 



3 In a Local Reference Frame centered on the detector D (see Fig. 1) the Nadir is the point on the 
celestial sphere which lies below the detector along its vertical axis, i.e., it is the point opposite to the 
detector Zenith. 



where [hi, h 2 ] is the Nadir angle interval of interest, 5 is the sampling function and 



T res = [ 2 dt5(hi < h{t) < h 2 ) (9) 

is the residence time, i.e., the time spent by the Sun in the Nadir angle interval. 

The values of < P e2 > and T res are functions of the Nadir angle interval, the time interval 
and of the sampling scheme described by the sampling function. Here it is assumed that 
the sampling i) is continuous over the time interval in which the Sun is in the relevant Nadir 
angle interval [h\, h 2 ], ii) is extended over an integer set of years, and iii) that the error in 
the reconstruction of the apparent Sun position at the time of the detection of each solar 
neutrino event is negligible. With these hypotheses 5(h\ < h < h 2 ) is 1 when h is inside 
the interval [h\, h 2 ] and zero otherwise. 

The choice of the Nadir angle interval is quite arbitrary, but a natural one is to compare 
data collected at Day (h > 90°) with data collected during the full Night (0° < h < 90°), 
and to separate neutrinos detected at night into Core neutrinos, 0° < h < h c , and Mantle 
neutrinos, h c < h < 90°, where h c = 33.17° is the apparent Nadir angle of the core/mantle 
boundary in the Stacey model. The three averages are labeled respectively: < P e2 > , 
< Pe2 > M and < P e2 > N , while the corresponding residence times are labeled T^ es , T^ s ; the 
Night residence time T^ es is identical to half an year for any detector. The day averaged 
probability < P e2 >d ay coincides with sin 2 6V- Finally, it is possible to consider also the 
Deep - Core average probability for which the Nadir angle interval in eq. (§) is h m ^ min < 
h < 2/3h c . The symbols Day, Night, Mantle, Core and Deep - Core will be used in what 
follows to indicate these particular solar neutrino event samples. 

In the absence of neutrino oscillations the number of solar neutrinos induced events 
in a given Nadir angle interval is proportional to T res for that interval. Table II lists a 
set of residence times for the Super - Kamiokande detector, showing that in one year, for 
approximately 14% of the total night time the Sun is seen behind the Earth core, so that 
about 7% of neutrinos cross the core. With a statistics of 10 4 events this corresponds to a 
relative statistical error of about 3.8%. 

The residence time is also relevant because it determines the weight of the contribution 
of each geophysical structure s, s =Night, Mantle, Core, Deep - Core, to < P e2 > N , the 
contribution being proportional to the ratio T^ esid /T y /2, where T y = 365.24 days. Thus, 
for the Super - Kamiokande detector the Core contribution in < P e2 > N has to be weighted 
by the factor 0.07 and it is no longer a dominant contribution even when < P e2 > c is much 
larger than < P e2 > M . 

Many types of systematic errors can affect the predictions for the D-N asymmetry. 
Among them the detector specifications such as detection efficiency, energy resolution, etc. 
are undoubtedly very important. They are neglected in the present study because they are 
not yet exactly known for the Super - Kamiokande detector, and moreover they may change 
during the various phases of data taking. In this case the only experimental feature to be 
taken into account is the energy threshold, while the most relevant error sources which have 
to be discussed are the numerical errors and the Sun's apparent motion reconstruction. Let 
us add that the error introduced by uncertainties in the detection cross section for the 



z/ e — e elastic scattering process is less than 1% [28] and is neglected. The same is true 
also for the detector location which is known with an accuracy better than 1 arcmin or 1 
over 10 4 M. 



The choice of a correct approximation for the apparent solar motion is a crucial point 
in order to reach the required accuracy, given the Earth model used in the calculations. 
The motion is quite complex and approximations are required to obtain an acceptable 
computation time but at the price of a certain error. There are two kinds of approximation 
made which are relevant for the Earth effect. In the first the Earth motion around the 
Sun is circular (COA), in the second elliptical (EOA). Neither the first nor the second are 
exact models, but EOA is very accurate for Earth effect prediction purposes. The analysis 
(see Appendix ^) shows that most of the errors introduced by COA and EOA are periodic 
in nature, so that in the turn of one year they should averaged out. However, given that 
P e2 is a function of the solar position and that the use of different sampling schemes (as 
Core, Mantle, Night, etc.) are equivalent to averaging over specific fractions of the time 
of the year only, it is not possible to test such approximations without accurate numerical 
experiments. These experiments show that indeed the error introduced by the simpler COA 
model can be close to (but does not exceed) 1%, and that a greater accuracy is reached by 
EOA. In this way the use of EOA for all the computations presented in this paper and in 
30| assures an accuracy better than 1% for each kind of sampling. The COA can be used 



when a large amount of simulations are required, allowing a substantial reduction of the 
computational time but at the cost of an increase in the error in < P e2 >, which can reach 
0.8%. 

In the Elliptical Orbit Approximation the averaged probabilities and residence times 
defined by (gj) have to be modified to take into account the change in the Sun - Earth 
distance with time, which causes a periodical change in the solar neutrino flux and in the 
Earth orbital velocity. To take both effects into account it is not enough to replace the 
relation between h and the time of the year t, but also the form of (|8|) has to be adapted. Let 
R(t) be the time dependent Sun - Earth distance expressed in units of the mean distance 
(the Astronomical Unit Rq = 1.4966 x 10 8 km). Then each instantaneous sample has to be 
weighted by the factor 1/R 2 {t) leading to: 



^ d 6(k < h(t) < h 2 ) 

'resJ Tl R 2 (t) 



<Pe2> = ThL* K ^ J P0W)), 



(10) 



T 2 

1 rP c I (lb 



6{h < h{t) < h 2 ) 



t, R\t) 



As it is discussed in Appendix |], the effect of the factor 1/R is partially compensated by 
the change in the residence time due to the change of the Earth orbital velocity during the 
year. 

3 The Instantaneous P e 2 Probability 

At the basis of each time averaged probability < P e2 > is the instantaneous one, P e2 . Many 
features of < P e2 > are a direct consequence of those present in P e2 . It is natural to begin 
the discussion of our numerical results with those obtained for P e2 . 

Figure 2 shows an example of the probability P e2 for sin 2 29y = 0.01 as a function of 
the resonance density p^ which is connected to E and the neutrino transition parameters 
Am 2 and sin 2 29y through the resonance condition: 

8 



p R Y e = 6.57 x 10 6 (Am 2 /E u )cos29 v , (11) 

where pR is expressed in gr/cm 3 , E v in MeV, Am 2 in eV 2 , Y e is the electron fraction per 
nucleon. We have for the Earth Y e = 1/2 with a rather high accuracy. For a fixed cos26y 
equation flTT| ) justifies the use of pr as an independent variable instead of E / Am 2 . The 
probability P e2 as a function of pr was computed solving the two-neutrino propagation 
equation |JT], [32] using a 4 4/l order, adaptive step-size Runge - Kutta algorithm adapted 
from p3[ . Each plot in Fig. 2 refers to a different Nadir angle. There are two main peaks, 
C (at pr = (10 — 11) gr/cm 3 ) and M. (at p R = (A — 6) gr/cm 3 ), associated with resonant 
transitions in the core (C) and in the mantle (A4). This correspondence is based on the 
following observations: i) a geological structures is relevant for the neutrino transitions (and 
therefore for P e2 ) when neutrinos cross it, ii) the structure affects the neutrino propagation 
when its density is the resonance region for the crossing neutrinos, iii) when one of the 
previous conditions are not fulfilled P e2 is close to its vacuum value sin 2 8y. When sin 2 26y 
increases, for instance, the resonance width also increases and it is possible to have a rather 
large P e2 for values of pr which do not correspond to any density crossed by neutrinos 
along their trajectory. In addition, the spherical symmetry of the Earth implies that solar 
neutrinos will cross the resonance region twice. This can produce interference terms in 
P e 2 leading to oscillatory dependence of P e2 on h for a fixed pr [|12]. They can also be 
responsible for the presence of multiple peaks in P e2 as a function of pr for a fixed h, 
associated with the resonances in a given geological structure. An example of this possibility 
is the splitting of the core peak C in Fig. 2(a) in two peaks, C/ and Co, shown in Fig. 2(c). 
Indeed, it is reasonable to expect the interference terms in P e2 to be significant in certain 
cases because the neutrino oscillation length in matter for solar neutrinos which undergo 
resonant transitions in the Earth is of the order of the Earth radius. 

It is instructive to see how the P e2 dependence on pr changes with the increasing of 
sin 2 20y. For sin 2 29y < 0.1 the qualitative properties of this dependences do not change 
significantly: P e2 is basically "rescaled" in accordance with the change of sin 2 29y and 
(P e2 — sin 2 9y) goes to zero rather quickly outside the resonance region. For larger values 
of sin 2 IQy the behavior is more complicate because new peaks appear, becoming more and 
more prominent as sin 2 29y increases. The "scaling" of P e2 is not precise, different parts 
of the plot scale differently, as is well illustrated by the C and M. peak heights. With 
the rise of sin 2 29y C increases faster than A4 reaching a maximum value and begins to 
decrease while M. is still increasing. As a consequence, at values of sin 2 26y £0.13 most of 
the D-N asymmetry is generated by the MSW effect in the core. This is not true in the 
large mixing angle region. For practical purposes the most important lesson which can be 
drawn from the P e2 plots is that to produce accurate predictions for the Earth effect, P e2 
has to be computed for resonance densities pr (or equivalently, for values of E / Am 2 ) which 
are largely outside the density range of the Earth geological structures. 

4 Averaging the P e 2 Probability 

Figures 3.1 - 3.18 depict the results of numerical calculations of < P e2 > and P© for the 
Super - Kamiokande detector in the case of solar v e transitions into an active neutrino 0, 

4 See appendix o for further details. 



v e — ► ^(r)- Each figure is divided into four "frames" which are labeled (a), (b), (c) and (d). 
Figures 3.1a - 3.18a display < P e2 > for Night (short-dashed line), Mantle (dotted line), 
Core (solid line) and Deep - Core (long-dashed line) as a function of pR and are the subject 
of the discussion in this Section. Comparing < P e2 > with P e2 shown in Fig. 2, it is possible 
to identify the peak due to the resonance in the mantle, M, present also in P e2 . Similarly, 
there are two very well defined peaks (or a peak, C, with two maxima) in < P e2 >, Cj (at 
larger p R ) and Co (at smaller p R ), which correspond to the resonance in the core. 

The presence of wiggles is evident for p R greater than the core density, especially at 
large mixing angles. These wiggles are dumped and disappear as pr increases. A curious 
feature is that in most of the cases the wiggles come in triplets (see Fig. 3.14a). The 
Log-Log plots of the < P e2 > probability (not reported here) show that for pr outside the 
range of the Earth density < P e2 > decreases as follows: 



!°g P^r£ -l)=A + B\og P R + ef(p R ), (12) 

where A, B and e are real constants and f{pR) is a dumped oscillating function which 
represents the wiggles. Since e is small (<10~ 2 ), /(/?#) can be neglected, while A and B can 
be determined by numerical means for any < P e2 >. Equation (|12|) is used to extrapolate 
< P e2 > at values of pr far from the Earth density. 

As Figs. 3. la - 3.18a illustrate, the separation in Core and Mantle samples is very effective 
in enhancing < P e2 >. The enhancement is particularly large for small mixing angles for 
which < P e2 > c can be up to six times bigger than < P e2 > M . This suggests the possibility 
of a corresponding enhancement of the e~-spectrum distortions as well as of the energy 
integrated event rate. The latter can make feasible the detection of the Earth effect even at 
small mixing angles. Further enhancement of the probability < P e2 > can take place if the 
size of the core bin is decreased, as shown by the Deep - Core averaged probability plotted 
with a long-dashed line. However, the reduction in the bin size reduces the residence time 
in the Deep - Core bin to less than 4% per year and it is not clear whether the probability 
enhancement can compensates the corresponding reduction in statistics. 

5 The Solar v e Survival Probability 

The connection between the peak structures in the probabilities < P e2 > and P 9 considered 
as functions of E/Am 2 , is determined basically by eq. ([5]). The presence of an enhancement 
in the averaged probability < P e2 > is not enough to ensure an enhancement in the D-N 
asymmetry. This is especially true at small mixing angles at which the Earth effect may 
produce reduction instead of increase of the v e flux. The latter happens when P & > 0.5, so 
that most of neutrinos coming from the Sun are v\ which are converted into v^, v T in the 
Earth. 

The upper parts ("windows") of the "frames" (b), (c) and (d) shown in Figs. 3.1 - 3.18 
are combined plots of P Q (dotted line), P e2 (dashed line) and P e (solid line) for Night (b), 
Core (c) and Mantle (d) as functions of E/Am 2 . The lower "windows" of each "frame" are 
plots of the ratio: 
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A P ^2^-& , (13) 

which mimics the asymmetry used by the Kamiokande collaboration in their discussion of 
the D-N effect. This quantity will also be called asymmetry. In some cases < P e 2 > is too 
small to be visible in the full scale plot, and so it is rescaled by a factor of 10. This is 
indicated by the label "< P e 2 > x 10" in the upper right corner of the corresponding figure. 
The purpose of these figures is to illustrate how each probability contributes in P and 
in the asymmetry Ap. One of the important features that has to be taken into account 
in the discussion of P is the position of the peaks in < P e 2 > with respect to the value 
of E/ Am? at which P = 0.5 because it controls the sign and the magnitude of the D-N 
effect. It follows from eq. (||) that the asymmetry is zero each time P = 0.5 irrespective 
of the value of < P e 2 >■ If < P e 2 > > sin 2 6V, but P > 0.5, the asymmetry is negative 
which means the Earth effect reduces the night event rate instead of increasing it. The 
positive Earth effects occurs when P < 0.5, i.e., when the peak of < P e 2 > falls inside 
the lower part of the "pit" of the MSW probability P . Figures 3.1 - 3.18 illustrate all 
these possibilities together with less relevant cases corresponding to < P e2 > = sin 2 9y and 

< P e 2 > < sin 2 #y- The presence of a zero Earth effect line together with cases in which 
the asymmetry Ap as a function of E/Am 2 , can be both negative and positive, can reduce 
the magnitude of the Earth effect. This reduction is especially important at small mixing 
angles. 

The increase of sin 2 29y within the small mixing angle solution region (3) does not change 
considerably the minimum and maximum values of P , but it affects significantly the width 
of the "pit" of Pq. At the same time the peaks in < P e2 > change mainly their height and 
width but not their position with respect to the E/Am 2 axis. For sin 2 2#y < 0.002 most 
of the relevant part of < P e 2 > is confined in the P > 0.5 region and the Earth effect 
is negative. With the increase of sin 2 29y the P "pit" becomes wider. Correspondingly, 
larger parts of the < P e 2 > peak enter the P < 0.5 region: depending on the value of 
E/Am 2 , the asymmetry Ap (the Earth effect) can be negative or positive. This continues 
until sin 2 29y = 0.006, for which most of the < P e 2 > peak is in the P < 0.5 region and 
the effect is mostly positive. This implies that even if < P e 2 > is relatively large, there will 
be a region at small mixing angles where the asymmetry is zero and the Core enhancement 
is not effective. 

The sequence of Ap plots allows to identify various peak components which we will 
denote by A, Ti + , HT , A being the left one (located at smaller values of E/Am 2 ), while 
7i + (H~) is located at larger E/Am 2 in the region where Ap > (Ap < 0). The A 
structure is an artifact of the presence of the adiabatic minimum in the probability P , 
minP & = sin 2 9y, which can lead to a relatively large asymmetry Ap even if < P e 2 > is 
small. The smallness of < P e 2 > in the A region means that that this region does not 
contribute significantly to the total event rate asymmetry. 

The bulk of the D-N asymmetry comes from 7i~ and 7i + structures, which are directly 
connected to the C and M peaks in < P e2 >. For sin 2 29y < 0.002 (Figs. 3.1 - 3.3) most 
of the peaks in < P e 2 > lie in the P > 0.5 region. Correspondingly, the asymmetry Ap is 
negative (only the peak 7i~ is present). It increases in absolute value as sin 2 29y increases 
and it is sensitive to the Core enhancement. As sin 2 29y grows further, the bulk of the 

< P e 2 > resonance part enters the region where P = 0.5, producing a suppression of the 

11 



asymmetry and of the Core enhancement, as is illustrated for sin 2 26V = 0.004 in Fig. 3.4. 
It is interesting to note that this crossing occurs at different sin 2 2Qy for the Night, Core 
and Mantle samples. As is shown in Fig. 3.4(c), the peak C for the Core is nearly cut in 
two parts by the P & = 0.5 line. This leads to the presence of both positive and negative 
Earth effect in the asymmetry Ap, while in the cases of Night and Mantle (Figs. 3.4(b) 
and 3.4(d)) the asymmetry is negative. At sin 2 29y = 0.006 (Fig. 3.5) the peak C is in the 
P Q < 0.5 region while the bulk of the effect of the Mantle is still in the P & > 0.5 region. 
The Core asymmetry is positive and it has a maximum value of 16.4%. For sin 2 29y > 0.008 
(Figs. 3.6 - 3.18) the effect is positive. All these details are influenced by the position of the 
Pq = 0.5 line with respect to the region of the < P e 2 > maxima. Given the fact that the 
latter is primarily a function of the Earth model, it cannot be excluded that certain aspects 
of the behavior described above may change somewhat with the change of the Earth model 
used in the calculations. 

Finally, the sequences of figures 3 show that the resonance in the core is the dominant 
source of the asymmetry in the region of the nonadiabatic solution, sin 2 29y < 0.013 while 
for sin 2 26y > 0.3 the bulk of the asymmetry is generated by the resonance in the mantle 
(the M peak) and the Core enhancement is only of the order of 30 -j- 40%. 

6 Conclusions 

In the present article and in ref . |3(| we have performed a detailed study of the D-N effect for 
the Super - Kamiokande detector assuming the MSW solution of the solar neutrino problem 
corresponding to solar v e transitions into an active neutrino, v e — > z>V r ). Among our aims 
were: i) calculation of the Earth effect with a sufficiently high precision which can match 
the precision of the data on the effect to be provided by the Super - Kamiokande detector, 
ii) calculation of the magnitude of the effect for a sufficiently large and representative set of 
values of the parameters from the "conservative" regions (^) and (^) of the MSW solution, 
iii) calculations of the one-year averaged D-N asymmetry for solar neutrinos crossing the 
Earth mantle, the core, the inner 2/3 of the core, and mantle + core (full Earth) for the 
chosen large set of values of Am 2 and sin 2 26y with the purpose of illustrating the magnitude 
of the enhancement of the effect which can be achieved by an appropriate selection of the 
data sample, iv) calculation of the one-year average deformations of the e~ spectrum by the 
Earth effect in the mantle, the core, and in the mantle + core for the selected representative 
set of values of Am 2 and sin 2 26V- In the present article we have used the elliptical orbit 
approximation for the Earth motion around the Sun to calculate the one year average solar 
v e survival probability and the D-N asymmetry in the probability for neutrinos crossing the 
Earth mantle, the core, the inner 2/3 of the core, and the mantle + core. The indicated 
sampling which can be done with the Super - Kamiokande detector, was considered in order 
to investigate quantitatively the possibility of enhancement of the D-N effect. We have 
found, in particular, that such an enhancement can be especially large at small vacuum 
mixing angles, sin 2 26y <0.013: the D-N asymmetry in the one-year averaged v e survival 
probability for neutrinos crossing the Earth core only can be a factor of up to six bigger 
than the asymmetry in the analogous probability for neutrinos crossing the core + mantle 



[Night data sample). This result persists in the corresponding event rate samples |30 ], thus 
suggesting that it may be possible to test the sin 2 26 v < 0.01 region of the MSW solution 
of the solar neutrino problem by performing selective D-N asymmetry measurements with 
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the Super - Kamiokande detector. 

We have studied also the accuracy of the circular orbit approximation (COA) utilized 
in many of the previous studies of the D-N effect, and compared it with the accuracy 
of the elliptical orbit approximation (EOA) we have used in the present analysis. The 
main conclusion of this study is that the largest difference in the results for the probability 
< P e2 > obtained using the COA and EOA is 0.8%. Given the fact that the EOA was chosen 
for the computations in this work, and that the EOA program actually includes most of the 
secondary celestial mechanics effects, the accuracy of the predictions for < P e2 > in what 
concerns the motion of the Sun is better (and probably much better) than 1%. 
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Appendix: The Solar Motion Approximations 

The calculation of the probability < P e2 > with a sufficient accuracy requires a relatively 
good approximation for the Sun's apparent motion which to leading order can be considered 
as a superposition of two strictly periodical motions, the orbital revolution and the Earth 
rotation, plus a set of small perturbations. The exact solution of such perturbation problem 
is complex and computationally demanding and therefore one has to select only those effects 
which are relevant for the chosen accuracy threshold e (ref. p7| |34| , |35| , |36|1). The various 



errors introduced by inadequate approximations are systematic in nature. They combine 
algebraically leading to subtle cancelation or enhancement effects which are difficult to 
propagate to the probability < P e2 > of interest. However, a practical approach to cut off 
all corrections which are not relevant is to set a threshold in the Sun's motion tracking 
accuracy of the order of some arcmin. This threshold may be justified by the fact that the 
Sun is not a point like neutrino source and the neutrino production region has an apparent 
angular diameter of about 3 arcmin (approximately one tenth of the solar diameter). The 
correctness of the indicated choice is confirmed by a more detailed analysis. 

Using the accuracy threshold it is possible to discriminate between the various assump- 
tions which can be used to simplify the apparent solar motion. As a result of such an 
analysis it becomes evident that two orbital approximations are sufficiently accurate for 
our purposes: the circular orbit approximation (COA) and the Elliptical Orbit Approxi- 
mation (EOA). In the COA the Sun's orbit is circular, while in the EOA the true elliptical 
Sun's orbit described by Kepler's laws is considered. There are two important effects which 
are neglected by COA and which instead are included in EOA: i) the time dependence 
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of the Sun - Earth distance, R(t), during the year, which leads to a change in the solar 
neutrino flux at the Earth, <J>, according to the 1/R 2 (t) law, ii) the time dependence of the 
Earth orbital velocity, v(t)<$, during the year due to the angular momentum conservation. 
The velocity f (£)© also changes according to the 1/R 2 (t) law. This implies a change in 
the residence time of the Earth, T res ffi , in a given orbital location, which is proportional to 
R 2 (t). Each of these two effects can introduce a systematic error of up to 3%. When taken 
together into account they compensate partially each other because the total neutrino rate 
for a given position of the Earth with respect to the Sun is proportional to the product 
T res ffi $ and therefore the total effect has to be rather small for a one year average Night 
sample. However, when other samples as Core or Deep - Core are considered cancelation 
at the 1% level cannot be guaranteed by pure heuristic arguments since even with an one 
year data taking the Core and Deep - Core sampling are equivalent to averaging over a 
fraction of the year (winter). In addition, for detectors located at positive latitudes over 
the tropical band, the Core sample is detected at winter when $ and v(t)$ are maximal. 

For the aforementioned reasons a set of numerical simulations were performed using 
both COA and EOA approximations and comparing the results for both < P e2 > and T res 
obtained with them. As already stated in the text, the main conclusion of this study is 
that the largest difference in the results for < P e 2 > obtained utilizing the COA and EOA is 
0.8%, and that the 3 arcmin error threshold provides a sufficient accuracy in the description 
of the solar motion. Let us note also that the program which make use of the elliptical orbit 
approximation takes into account many other celestial mechanics secondary effects as well, 
so that its accuracy is better than few tenth of an arcmin. Since all the results presented 
in this paper are obtained with this program, it is evident that the accuracy of < P e 2 > 
predictions here presented is far better than 1%. 

A number of tests where also performed to assure that the accuracy of the numerical 
calculations is sufficiently high. For instance, calculated solar positions were compared with 



positions given in the 1996 edition of the Nautical Almanac p6| . These tests showed that 
the numerical errors are not greater than 10 _3 % and therefore are negligible. 
Finally, let us note that eq. @ can also be written using the ratio 

i Q-* res 



T res (h<90°) dh 

i.e., the probability to receive a neutrino from the Sun (of any flavour or energy) when 
it is inside the Nadir angle interval h,h + dh. In this case the expression for the average 
transition probability for the Nadir angle interval [h\, hW becomes: 

Although it looks more attractive theoretically, this formulation is unpractical due to the 
presence of 365 singularity points in the dT res /dh function, as is illustrated in Fig. 4. These 
singularities are due to the fact that when the Sun's Nadir angle reaches its minimum value, 
dh(t)/dt vanishes while dT res /dh ~ l/(dh(t)/dt) diverges. Obviously, these singularities are 
integrable but they make rather complicated the use of eq. (A.l) for sufficiently accurate 
calculations. 
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Tables 



Table I 


List of the A r 


n 2 and sin 2 26 


sample values 






Am 2 (eV 2 ) 




sin 2 26 v 


I II 


III 


IV 


0.0008 


9E-6 7E-6 


5E-6 




0.0010 


9E-5 7E-6 


5E-6 




0.0020 


1E-5 7E-6 


5E-6 




0.0040 


1E-5 7E-6 


5E-6 




0.0060 


1E-5 7E-6 


5E-6 




0.0080 


1E-5 7E-6 


5E-6 




0.0100 


7E-6 5E-6 






0.0130 


5E-6 






0.3000 


1.5E-5 2.0E-5 


3.0E-5 


4.0E-5 


0.4800 


3E-5 5E-5 






0.5000 


2E-5 






0.560 


1E-5 






0.600 


8E-5 






0.700 


3E-5 5e-5 






0.770 


2E-5 






0.800 


1.3E-4 






0.900 









Table II. Residence times for the Super - Kamiokande detector. 





Nadir Angle h 








Sample s 


from to 


Jp x 100 
1 y 


ipr xioo 


4^x100 

K-E 


Deep - Core 


0° 25.98° 


3.94 


7.86 


5.0 


Core 


0° 33.17° 


7.11 


14.18 


3.8 


Mantle 


33.17° 90° 


43.03 


85.82 


1.5 


Night 


0° 90° 


50.14 


- 


1.4 


Note: 




T y = 365.243 days. 




The integration is extended o^ 


rer one year from day to day 365 


o-(7Ze) is the 1 year mean evei 


it rate statistical uncertainty for TZe,o — 10 4 i^/year 
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Figure Captions 

Figure 1. The geometry of the Earth effect: D and denote the detector and the Sun, 
h is the Sun's Nadir angle, P® is the Earth radius (6371 km) and R cor e is the core radius 
(3486 km). The small arrows point in the directions of the detector Zenith and Nadir. 



Figure 2. The dependence of the probability P e 2 on pn for sin 2 2#y = 0.01. The five 
plots are obtained for O.lgr/cm < pn < 30.0 gr/cm and five different solar neutrino 
trajectories in the Earth: a) center crossing (h = 0°), b) winter solstice for the Super - 
Kamiokande detector (h = 13°), c) half core for the Super - Kamiokande detector (h = 23°), 
d) core/mantle boundary (h = 33°), e) half mantle (h = 51°). 



Figure 3. The probabilities < P e 2 > and P® as functions of E/Am 2 . Each figure represents 
P®, < P e 2 >, P® an d Ap for the value of sin 2 29y indicated in the figure. The frames (a) 
show the probability < P e 2 > as a function of the resonance density pn for the Night (short 
- dashed Line), Mantle (dotted line), Core (full line) and Deep - Core (long - dashed line) 
samples. The frames (b), (c) and (d) represent the probabilities P® (Dotted Line), < P e 2 > 
(Dashed Line), P® (Full Line) corresponding to the Night, Core and Mantle samples as 
functions of E u /Am 2 . In some cases < P e 2 > is multiplied by a factor of 10, which is 
indicated in the corresponding figures. 



Figure 4. Distribution of the differential residence time dT res /dh° for the Super - 
Kamiokande detector. The Nadir angle is expressed in degrees. The upper figure depicts 
the distribution for h < 89°, while the lower figure shows an enlargement of the distribution 
for 40° < h < 44° . The lower figure shows also that the small wiggles in the upper figure are 
due to a set of singularities not well sampled by the computing algorithm. The peaks in the 
lower figure have a finite height only because of limitations in the numerical computation 
accuracy. 
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